Combinatorics and Boson normal ordering: A gentle introduction 
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We discuss a general combinatorial framework for operator ordering problems by applying it to the 
normal ordering of the powers and exponential of the boson number operator. The solution of the 
problem is given in terms of Bell and Stirling numbers enumerating partitions of a set. This frame- 
work reveals several inherent relations between ordering problems and combinatorial objects, and 
displays the analytical background to Wick's theorem. The methodology can be straightforwardly 
generalized from the simple example given herein to a wide class of operators. 



I. INTRODUCTION 

Hilbert space constitutes the arena where quantum 
phenomena can be described. One common realization is 
Fock space, which is generated by the set of orthonormal 
vectors \n) representing states with a specified numbers of 
particles or objects. A particular role in this description 
is played by the creation and annihilation a operators 
representing the process of increasing and decreasing the 
number of particles in a system, respectively. We con- 
sider operators that satisfy the boson commutation rela- 
tion [a, a^} = 1 describing objects obeying Bose-Einstein 
statistics, for example, photons or phonons. The fact 
that the operators a and do not commute is proba- 
bly the most prominent characteristic of quantum the- 
ory, and makes it so strange and successful at the same 



time 
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In this paper we are concerned with the order- 
ing problem which is one of the consequences of non- 
commutativity. This problem derives from the fact that 
the order in which the operators occur is relevant, for ex- 
ample, a^a 7^ aa^ = a^a+ 1. The ordering problem plays 
an important role in the construction of quantum me- 
chanical operators. The physical properties of differently 
ordered operators may be quite distinct, which we can 
see by considering their expectation values. An analysis 
of operator matrix elements reveals their physical prop- 
erties observed as probabilities. There are two sets of 
states of primary interest in this context: number states 
\n) and coherent states \z). The latter, defined as eigen- 
states of the annihilation operator a, play an important 
role in quantum optics^£&£ and in the phase space for- 
mulation of quantum mechanics^ 

The calculation of the number or coherent state expec- 
tation values reduces to transforming the original expres- 
sion to the normally ordered form in which all annihila- 
tion operators are to the right. In this form the evalua- 
tion of the matrix elements is immediate. The procedure 



is called normal ordering 



Although the process is 



clear and straightforward, in practice it can be tedious 
and cumbersome when the expression is complicated, and 
is even less tractable when we consider operators defined 
by an infinite series expansion. It is thus desirable to 
find manageable formulas or guiding principles that lead 
to solutions of normal ordering problems. 

In this paper we present a general framework that is 
applicable to a broad class of ordering problems. It ex- 
ploits the fact that the coefficients emerging in the nor- 
mal ordering procedure appear to be natural numbers 
which have their origin in combinatorial analysis. In the 
simplest case of powers or the exponential of the num- 
ber operator N = a^a these are Stirling and Bell numbers 
which enumerate partitions of a set £ We use this example 
to illustrate a systematic approach to the ordering prob- 
lem. The general methodology is to identify the problem 
with combinatorial structures and then resolve it using 
this identification. The solution may be found with the 
help of the Dobinski relation, 10,11 which is a very effective 
tool and is straightforwardly applicable to a wide range 
of ordering problems. 

As a byproduct of this methodology we obtain a sur- 
prising relation between combinatorial structures and op- 
erator ordering procedures. This relation is especially in- 
teresting because the objects involved in the problem can 
have clear combinatorial interpretations (for example, as 
partitions of a set). The expectation is that this remark- 
able interrelation will shed light on the ordering problem 
and clarify the meaning of the associated abstract oper- 
ator expressions. 

The framework we present is an example of a fertile in- 
terplay between algebra and combinatorics in the context 
of quantum mechanics. It employs only undergraduate 
algebra and is not as yet a standard feature of quantum 
mechanics textbooks. 

The paper is organized as follows. Section [IT] briefly re- 
calls the concept of Fock space and introduces the normal 
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ordering problem. The main part containing the connec- 
tion to combinatorics is given in Sec. IIII1 It illustrates 
the methodology by discussing in detail the solution of 
a generic example. Some applications are provided in 
Sec. IIV1 In Sec. [V] we point out extensions of this ap- 
proach and suggest further reading. 



These states are normalized, (z\z) = 1, but are not or- 
thogonal and constitute an overcomplete basis in the 
Hilbert space.— Coherent states have many useful prop- 
erties which are exploited in quantum optics and in other 
areas of physics £d£JL2JL 



II. OCCUPATION NUMBER 
REPRESENTATION 

A. States and operators 

We consider a pair of one mode boson annihilation 
a and creation operators satisfying the conventional 
boson commutation relation 



[a, a 1 ] = 1. 



(1) 



The operators a, a^, and 1 generate the Heisenberg alge- 
bra. 

The occupation number representation arises from the 
interpretation of a and as operators annihilating and 
creating a particle in a system. From this point of view 
the Hilbert space TL of states is generated by the number 
states |n), where n = 0,1,2,... counts the number of 
particles, or objects in general. The states are assumed 
to be orthonormal, (m\n) = 5 m . n , and constitute a basis 
in H. This representation is usually called Fock space. 

The operators a and satisfying Eq. fl} may be re- 
alized in Fock space as 



a \n) = \fn\n — 1) 



*t| n ) = Vn+l|n+l). (2) 



The number operator N, which counts the number of 
particles in a system, is defined by 



N\n) = n\n) 



(3) 



and is represented as N = a^a. It satisfies the commuta- 
tion relations 



[a,N] 



,t 



N] 



-,t 



(4) 



The algebra defined by Eqs. ((T|) and ((4]) describes objects 
obeying Bose-Einstein statistics, for example, photons or 
phonons. It is sometimes called the Heisenberg- Weyl al- 
gebra, and occupies a prominent role in quantum optics, 
condensed matter physics, and quantum field theory. 

The second set of states of interest in Fock space are 
the coherent states \z). They are defined as the eigen- 
states of the annihilation operator 



B. Normal ordering: Introduction 

The noncommutativity of the creation and annihilation 
operators causes serious ambiguities in defining operator 
functions in quantum mechanics. To solve this problem 
the order of the operators has to be fixed. An important 
practical example of operator ordering is the normally 
ordered form in which all annihilation operators a stand 
to the right of the creation operators a) . We now define 
two procedures on boson expressions yielding a normally 
ordered form, namely, normal ordering and the double 
dot operation.— 

By the normal ordering of a general expression F(a i , a) 
we mean F^(a\ a) which is obtained by moving all the 
annihilation operators a to the right using the commuta- 
tion relation of Eq. ([1]) . This procedure yields an opera- 
tor whose action is equivalent to the original one, that is, 
F^(a\a) = F(a^ , a) as operators, although the form of 
the expressions in terms of a and may be completely 
different. 

The double dot operation :F(a',a): consists of apply- 
ing the same ordering procedure but without taking into 
account the commutation relation of Eq. |T]), that is, 
moving all annihilation operators a to the right as if they 
commuted with the creation operators . This nota- 
tion, although widely used, is not universal.— We observe 
that in general this procedure yields a different operator 
F{a\a)^:F{a\a): ^ 

In addition to the fact that these two procedures yield 
different results (except for operators that are already in 
normally ordered form), there is also a practical differ- 
ence in their use. That is, although the application of 
the double dot operation is almost immediate, for the 
normal ordering procedure some algebraic manipulation 
of the non-commuting operators a and a' is needed. Here 
is an example of both procedures in action: 



I 4- normal ordering 

aa'aaa'a — 

[a,at] = l 



(a t ) 2 a 4 + 4a t a 3 + 2 



to the left, a to the right 



aa^aaa^a 



double dot 



a. a 1 commute 
(like numbers) 



a) a) aaaa . 



a\z) = z\z), (5) 

where z is a complex number (the dual relation is (z\a< = 
z*{z\). These states take the explicit form 



e 2 



\4 



E 



(6) 



In general we say that the normal ordering problem for 
F(a^, a) is solved if we can find an operator G(a^ , a) for 
which the following equality is satisfied 



F{a\a) = :G(a\a): 



(7) 



The normally ordered form has the merit of enabling im- 
mediate calculation of an operator's coherent state ele- 
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ments which reduce, by virtue of Eq. ([5]), to substituting 
a — > z and o) — > z* in its functional representation, that 
is, 



(z|:G(at,o):|z) =G(z*,«). 



(8) 



Thus, by solving the normal ordering problem of Eq. ([7]), 
we readily obtain 

{z\F(a\a)\z) = G(z*,z) . (9) 

This procedure may be illustrated in the example 

{z\aa) aaa) a \z) = (z|(a t ) 2 a 4 + 4a t a 3 + 2a 2 \z) 
= {z*) 2 z 4 + Az*z 3 + 2z 2 . 

In brief, we have shown that the calculation of coher- 
ent state matrix elements reduces to solving the normal 
ordering problem. The converse statement is also true; 
that is, if we know the coherent state expectation value 
of the operator, say Eq. ((9]), than the normally ordered 
form of the operator is given by Eq. 

A standard approach to the normal ordering problem 
is to use Wick's theorem.— In our context, this theorem 
expresses the normal ordering of an operator by applying 
the double dot operation to the sum of all possible ex- 
pressions obtained by removing pairs of annihilation and 
creation operators where a precedes a\ called contrac- 
tions in analogy to quantum field theory, for example 



aa^aaa/a = 



aoJaaoJa : 



no pair removed 

v aaa' a + daJaaeP a + aa'dad'a + aa/add'a : 



1 pair removed 

tfda$a + ^d ] a^a : = (a r ) 2 a i + 4 a)a a + 2 a 2 



2 pairs removed 

This procedure may involve a large number of steps. For 
polynomial expressions this difficulty may be overcome 
by using computer algebra, although this use does not 
provide an analytic structure. For nontrivial functions, 
such as those having infinite expansions, the problem still 
remains open. 

One approach to the problem relies on the disentan- 
gling properties of Lie algebraic operators and applica- 
tion of the Baker-Campbell-Hausdorff formula. Here is a 
standard example: 



e A(a+a f ) _ g A 2 /2 . e A(a+a f ). 



(10) 



However, the use of this kind of disentangling property 
of the exponential operators is restricted in practice to 
quadratic expressions in boson operators^ 

Another method exploits the recurrence relations and 
solves the normal ordering problem by use of combina- 
torial identities. 9,17 This promising approach was the in- 
spiration for the systematic combinatorial methodology 
which is presented in this article. 



III. GENERIC EXAMPLE: STIRLING AND 
BELL NUMBERS 

A. Normal ordering: Combinatorial setting 

We consider the number operator N — a^a and seek 
the normally ordered form of its nth power [n > 1). We 
write the latter as 



(o)a) n = J2S(n,k) (a t ) fe a fc , 



(11) 



fc=i 



which uniquely defines the integer sequences S(n, k) 
for k = 1 . . . n; these sequences are called the Stirling 
numbers i 11 ^ Information about this sequence for each n 
may be captured in the Bell polynomials 



B(n, x) = Y,S(n,k)x k . (12) 
fc=i 

We also define the Bell numbers B(n) = B(n, 1) as 

n 

B(n) = J2S(n,k). (13) 



k=l 



Instead of operators a and we may equally well insert 
into Eq. (fTl"]) the representation of Eq. (fT]) given by the 
operator X defined as multiplication by x, and by the 
derivative D 



19 



dx 



x, 



D. 



(14) 



This substitution does not affect the commutator of 
Eq. (fT]), that is, [D, X] — 1, which is the only property 
relevant for the construction. Therefore in this represen- 
tation Eq. (jlip takes the form 



(ID)" = ^S{n,k)X k D k . 



(15) 



k=l 



B. Combinatorial analysis 

In the following we discuss the properties of the Stir- 
ling and Bell numbers^ For that purpose we use ele- 
mentary methods of combinatorial analysis based on a 
versatile tool known as the Dobihski relation ! 10 ! 11 The 
latter is obtained by acting with Eq. (fTS"]) on the expo- 
nential function e x = J^kLo T\ yielding 



5> 

k=0 



kl 



■^2S(n,k): 



(16) 



fc=i 



[xx all eqs. must be numbered xx] If we recall the defini- 
tion of the Bell polynomials Eq. (fT2|) . we obtain 



(17) 



k n 

B(n,x) =e~ x ^2^x k , 

fc=0 
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which is the celebrated Dobinski relation ^ ! 11 It is usually 
expressed in terms of the Bell numbers, which are given 
by 



B(n) 



k=0 



k" 



(18) 



We observe that both series are convergent and express 
the integers B(n) or polynomials B(n,x) in a nontriv- 
ial way. To mention one of the many applications notice 
that k n in Eq. (fTT)) may be replaced by the integral repre- 
sentation k n = J Q dXX n S(X — k). If we change the order 
of the sum and integral (allowable because both are con- 
vergent), we obtain the solution to the Stieltjes moment 
problem for the sequence of Bell polynomials 



where 



B(n, x) 



W X (X) 



dX W x (X)X ri 



gj(A 

k=0 



k\ 



—x k 



(19) 



(20) 



is a positive weight function located at integer points 
and is called a Dirac comb.— Note that Eq. (|20|) may be 
identified with the Poisson distribution with mean value 
equal to x. 

A very elegant and efficient way of storing and tackling 
information about sequences is attained through their 
generating functions. — The exponential generating func- 
tion of the polynomials B(n,x) is defined as 



G(X,x) = Y i B{n,x) 



A™ 



(21) 



It contains all the information about the Bell polynomi- 
als. If we use the Dobinski relation, we may calculate 
G(X, x) explicitly. We substitute Eq. (H} into Eq. (|2Tjl . 
change the summation order,— and then identify the ex- 
pansions of the exponential functions to obtain 



oo oo 



G(X,x) 



■ \ " \ h n — — 

n=0 k=0 
00 ^.k cv \n 

E — k n — 
k\ ^ n\ 

k=Q n=0 



E 

fc=0 



k\ 



(22a) 
(22b) 
(22c) 



Thus, the exponential generating function G(A, x) takes 
the compact form 



G(A, x) 



Note that in the context of the weight function of Eq. (j2"D|) 
G(A, x) is the moment generating function of the Poisson 
distribution with the parameter x. 



An explicit expression for the Stirling numbers S(n, k) 
may be extracted from the Dobinski relation. Note that 
in Eq. (TTTj) the relevant series may be multiplied together 
using the Cauchy multiplication rule to yield 



oo i OO 

Bfcx) = 5>i)'££* w Tr 



1=0 
oo k 



k=0 



(24a) 



(24b) 



fe=0 j=l 

By comparing the expansion coefficients in Eq. (|24b[) with 
Eq. (fT2"|) , we obtain 



S(n, k) 



1 k 
kl E 



(-l) k - j 3 n , 



(25) 



which yields an expression for S(n, k). 

If we use any of the above standard formulas Eqs. 
and (fT3|) for the Stirling or Bell numbers, we can easily 
calculate them explicitly. We remark that many other in- 
teresting results may be derived by straightforward ma- 
nipulation of the Dobinski relation Eq. (fT7|) or the gener- 
ating function Eq. ((23|) . see Appendix IA"1 23 Some of these 
numbers are given in Table: 



S(n, k), 



Kk <n 



B(n) 



n 
n 
n 
n 
n 
n 
n 
n 



1 

1 1 

1 3 1 

1 7 6 

1 15 25 

1 31 



1 

10 

90 65 



1 

15 



1 63 301 350 140 21 1 
1 127 966 1701 1050 266 28 1 



1 

2 
5 

15 

52 
203 
877 
4140 



C. Normal ordering: Solution 

We return to normal ordering. By using the proper- 
ties of coherent states in Eqs. ([T])-©, we conclude from 
Eqs. (jlip and (|12p that the diagonal coherent state ma- 
trix elements generate the Bell polynomials 17 



(z\(aU) n \z) =B(n,\z\ 2 ). 



(26) 



If we expand the exponential e Aa a and take the diagonal 
coherent state matrix element, we find 



(z\e^ a \z) 



OO 

E< 

n=0 



z\{(Ja) n \z) 



X" 



E B (™'N 2 ) 



A" 



n=0 



(27) 

(23) We observe that the diagonal coherent state matrix ele- 
ments of e Aata yield the exponential generating function 
of the Bell polynomials (see Eqs. (JUl) and (f2"3"]0 



(z\e XaTa \z) 



(28) 



Equations (|26|) and (|28|) allow us to read off the normally 
ordered forms 



and 



(a f a)" = :B(n,a^a): , 

e Aa f a _ . e a t a(e A -l). 



(29) 
(30) 



S(3,3)=l 



Notice that the normal ordering of the exponential of 
the number operator a^a amounts to a rescaling of the 
parameter A — > e A — 1. We stress that this rescaling is 
characteristic for this specific case only and in general the 
functional representation may change significantly (see 
Sec.|V|. Just for illustration we give results that can be 
obtained by an analogous calculation 24 



e A < at ) 2a = :exp 



/A£aW\. 
^ 1-Aa+ J ■ ' 



s A(at) 2 a 2 



-atav^oo An(n-l) (a' a)" , 
Z^n=0 e n! ' 



(31) 
(32) 



A. ft 





U da] 



S(3,2)=3 

FIG. 1: Illustration of Stirling numbers S{n,k) enumerating 
partitions of a set of n = 3 distinguishable marbles (white, 
gray and black) into k = 1, 2, 3 su bsets. 



Equations (f2"9"| and (f3"0")) provide an explicit solution to 
the normal ordering problem for powers and the exponen- 
tial of the number operator. This solution was obtained 
by identifying the combinatorial objects and resolving 
the problem on that basis. Furthermore, this surprising 
connection opens a promising approach to the ordering 
problem through its combinatorial interpretation. 



D. Combinatorial interpretation: Bell and Stirling 
numbers 

We have defined and investigated the Stirling and 
Bell numbers as solutions to the normal ordering prob- 
lem. These numbers are well known in combinatoric o 1 1 i 18 
where the S(n, k) are called Stirling numbers of the sec- 
ond kind. Their original definition is given in terms of 
partitions of a set; that is, the Stirling numbers S(n, k) 
count the number of ways of putting n different objects 
into k identical containers (none left empty), see Fig. Q] 
The Bell numbers B(n) count the number of ways of 
putting n different objects into n identical containers 
(some may be left empty). From these definitions the re- 
currence relation of Eq. (|A1|) may be readily obtained and 
further investigated from a purely combinatorial view- 
point. This formal correspondence establishes a direct 
link to the normal ordering problem of the number oper- 
ator. As a result we obtain an interesting interpretation 
of the ordering procedure in terms of combinatorial ob- 
jects. We remark that other pictorial representations can 
be also given, for example, in terms of graphs 25 or rook 
numbers.— 

In conclusion we point out that this method may be re- 
versed; that is, certain combinatorial families of numbers 
may be given an algebraic interpretation in the quantum 
mechanical context. 



IV. SOME APPLICATIONS 

A. Quantum phase space 

A curious application of the coherent state represen- 
tation is found in the phase space picture of quantum 
mechanics. For the conjugate pair of observables q = 
(a) + a)/y/2 and p = i(a i ' — a)/V2, related to the posi- 
tion and momentum operators of a particle or to quadra- 
tures of the electromagnetic field, coherent state expecta- 
tion values have the simple form (z|g|z) + i(z\p\z) = \[2z 
and minimize the uncertainty relational In this sense the 
coherent state \z) for z = (q + ip)/y2 may be interpreted 
as the closest quantum approximation to the classical 
phase state (q,p). These properties are used to construct 
the quantum analog of phase space through the Husimi 
distribution, denoted by Q(q,p), which for the quantum 
state described by the density matrix p is defined as^^ 

Q(q,p) = ^(z\p\z). (33) 

Q(q,p) is interpreted as the probability density for the 
system to occupy a region in phase space of width 
Aq = Ap = \J\j1 centered at (q,p) which on experi- 
mental grounds refers to obtaining the result (q,p) from 
an optimal simultaneous measurement of q and p. Such 
measurements in quantum optics are obtained using the 
technique of heterodyne detection. 

This construction of the quantum phase space analog 
raises the problem of efficiently calculating the coherent 
state expectation values of an operator which, as we have 
seen in Sec. Ill Bt is in practice equivalent to its normal 
ordering. Hence ordering techniques are important for 
practical use. 

As an illustration we observe that from Eq. (|28|) we 
can readily derive the explicit expression for the Husimi 
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distribution of the quantum harmonic oscillator in ther- 
mal equilibrium. For the Hamiltonian H = a^a + 1/2 the 

density matrix p of a thermal state is e~ /3a ^ a /Z : where 
Z = Tre-^ a = 1/(1- e' 13 ) and = l/k B T. Thus from 
Eqs. ([33]) and (28|) we obtain 



Q(q,p) 



2tt 



(1 



^)e<- e ~ 1){q +p )/2 . 



(34) 



It is instructive to compare this quantum phase space 
distribution with its classical analog. The correspond- 
ing Hamiltonian for the classical harmonic oscillator is 
H c \ = {q 2 +P 2 )/2, and the probability distribution in 
the thermal state is P c \(q,p) = e - ^ 9 +p ^ 2 /Z c i, where 
Zd = / e -^« 2+ P 2 )/ 2 dqdp = 2tt/0. Finally, we obtain 



P cl (q,p) = J-/? e -/% 2 +/)/ 2 . 

Z7T 



(35) 



In both cases we obtain gaussians. However, observe that 
the quantum distribution of Eq. ([53]) is wider than its 
classical analog of Eq. (|35|) . It is explained by additional 
fluctuations due to the uncertainty relation that are in- 
herent in quantum mechanics. For — > 0, that is, for 
large temperatures, the quantum distribution of Eq. (|34|) 
correctly goes to the classical distribution in Eq. |35|) . 

An analogous analysis can be done for the whole spec- 
trum of models described by Hamiltonians constructed 
in the second quantization formalism, provided the nor- 
mally ordered form of the operators is known. Section IHT1 
discusses the methodology which is applicable to a wide 
set of problems, for example, the optical Kerr medium as 
in Eq. (|3"2"|) and the open system described by Eq. (|3"Tj) . 
See Sec. |V] for a discussion of the range of applicability. 



B. Beyond the Wick theorem 



As mentioned in Sec. IIIBI the standard approach to 
normal ordering through Wick's theorem reduces the 
problem to finding all possible contractions in the op- 
erator expression. In practice, the process may be te- 
dious and cumbersome to perform, especially when a 
large number of operators are involved. Hence system- 
atic methods, like the one described in Sec. IIII1 are of 
importance in actual applications. 

To complete the picture we will show how to connect 
Wick's approach to the combinatorial setting described 
in this paper. The bridge is readily provided by the in- 
terpretation of Stirling numbers as partitions of a set, as 
given in Sec. IIIIDI To see the link we consider a string 
a^a a^a . . . a^a consisting of n blocks a^a which we la- 
bel from 1 to n starting from the left, thus obtaining n 
distinguishable objects 



a'a a'a 



a'a 



@© ••■ ® 



Then each choice of contraction in Wick's theorem 
uniquely divides this set into classes such that objects 



in the same class are connected by contractions between 
their operator constituents. See the following examples 
for illustration 



A a)a a)a a'a A <-> © © © © © <-> @© 



t 1 't' 't t 1 't 



tr-jtit l t n t 

"i 'ft rt'n n < n n <ri n <i 



-Fn lr 



a'a a'a a'a a'a a'a ^ © © © @ © @@ 



@©@ 


®@ 




©@ 


©@© 



Observe that this construction may be reversed and thus 
provides a one-to-one correspondence between operator 
contractions in (a'a) n and partitions of the set of n 
objects^ In this way the contractions of Wick's theo- 
rem may be seen as partitions of a set, providing the link 
to the combinatorial framework presented in this paper. 
The methodology of Sec. [TH] offers an alternative perspec- 
tive on the normal ordering problem and unlike Wick's 
approach, exposes its analytical structure, thus yielding 
practical calculational tools. 



C. Operator identities 

Manipulations of differently ordered operators often 
lead to interesting operator identities. For example, tak- 
ing the limit A — > — oo in Eq. (|30p. yields an interesting 
representation of the vacuum projection operator^ 



|0)(0|=:e- 



(36) 



Equation (f3"6"]) leads to a coordinate representation of the 
squeezing transformation 5(A) = e^ A ° ~ Xa which 
is extensively used in quantum optics^ It may be ob- 
tained using the technique of integration within an or- 
dered product, 30 yielding 



5(A) = e 



- p -A/2 



dq\e- X q)(q\ 



(37) 



which offers an interpretation of 5(A) as an explicit 
squeezing of the quadrature. 



V. SUMMARY AND OUTLOOK 

We have presented a combinatorial framework for op- 
erator ordering problems by discussing a simple example 
of the powers and exponential of the number operator 
a'a. We have provided a general method for relating 
normally ordered operator expressions to combinatorial 
objects and then solved the problem from that view- 
point. The solution involved using the representation of 
the Heisenberg algebra in terms of operators on the space 
of polynomials and then applying the Dobihski relation, 
which provides the exponential generating function and 
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explicit expressions. This approach provides effective cal- 
culation^ tools and also exposes the analytic structure 
behind the Wick theorem. 

One advantage of this methodology is that it can be 
straightforwardly generalized to a wide class of opera- 
tor expressions.— The simplest examples are provided 
by the powers and exponentials of (a^Ya and (a^) r a s 
with r and s integers. 31 It may be further extended to 
investigate the normal ordering of boson monomials 25 in 
the form (a'Y n a Sn . . . (a <! ) r2 a S2 (a' !! Y 1 a 31 and more gen- 
erally homogeneous boson polynomials^ 2 - that is, linear 
combinations of boson expressions with the same excess 
of creation over annihilation operators^ Further devel- 
opment of the method applies to the ordering of general 
operators linear only in the annihilation (or creation) op- 
erator, that is, q(a^)a + v(aJ), where q(x) and v(x) are 
arbitrary functions. The exponential of an operator of 
this type constitutes a generalized shift operator and the 
solution is in the class of Sheffer polynomials.— In all of 
these cases the use of the Dobihski relation additionally 
provides a solution of the moment problem,— as well as a 
wealth of combinatorial identities for sequences involved 
in the result (including their deformations^). 

Ordering problems are naturally inherent in the alge- 
braic structure of quantum mechanics. It is remarkable 
that they may be described and investigated using ob- 
jects having a clear combinatorial interpretation. For the 
generic example considered here these are partitions of a 
set. For more complicated expressions the interpretation 
can be provided by introducing correlations between el- 
ements or by using a graph representation. 

APPENDIX A: COMBINATORIAL IDENTITIES 

We enumerate some properties of the Stirling and 
Bell numbers defined in Sec. IIII^ 5 - The reader is invited 
to check these relations by straightforward manipula- 
tion of the Dobihski relation or exponential generating 



function^ 

The recurrence relation for the Stirling numbers is 

S(n + l,k) = kS(n,k) + S(n,k-l), (Al) 

with the initial conditions S(n, 0) = <5„.o and S(n, k) = 
for k > n. 

The Bell polynomials may be shown to satisfy the fol- 
lowing recurrence relation 

B(n+l,x) = xf2( 1 f)B(k,x), (A2) 

k=0 ^ ' 

with B(0, x) = 1. Consequently for Bell numbers we have 

B(n+l)=ELo©#); 

The following exponential generating function of the 

Stirling numbers S(n, k) is sometimes used in applica- 
tions 

00 \ n (p x i\ k 

n=k 

In addition, the Stirling numbers S(n, k) may be inter- 
preted as the connection coefficients between two sets x n 
and x—, n = 1, 2, . . ., where x— = x ■ (x— 1) . . . (x — n+ 1) 
is the falling factorial; that is, they represent a change of 
basis in the space of polynomials 

n 

x n = S(n, k) xK (A4) 

k=l 

We also note that the Bell polynomials belong to the 
class of Sheffer polynomials^ which in particular share 
an interesting property called the Sheffer identity (note 
the resemblance to the binomial identity) 

B(n, x + y) = Q B ( k > v) B ( n - k,x). (A5) 
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